Photochemical Reaction of Ketoprofen with Proteinogenic Amino Acids.

Ketoprofen (KP) is one of the most popular nonsteroidal anti-inflammatory drugs; however, drug-induced photosensitivity of KP has been reported as a serious adverse effect. KP incorporated into a protein can produce an allergen under UV irradiation, which causes drug-induced photosensitivity. The photochemistry of KP with 20 kinds of proteinogenic amino acids in phosphate buffer solutions at pH 7.4 was studied by transient absorption spectroscopy. The KP carboxylate anion (KP-) gave rise to a carbanion via a decarboxylation within a laser pulse, and the carbanion yielded 3-ethylbenzophenone ketyl biradical (3-EBPH) through a proton transfer reaction. Twelve kinds of proteinogenic amino acids obviously accelerated the reaction. Structural information on the complexes of KP docked in the binding sites of human serum albumin (HSA) was obtained by molecular mechanics (MM) and molecular dynamics (MD) calculations. The photochemical reaction of KP- with amino acid residues in HSA was discussed on the basis of the experimental and calculational results. The information on the reactivity of KP with the amino acids and the stable structures of the KP-HSA complexes should be essential for understanding of the initial step for drug-induced photosensitivity.

Thermal and Nonthermal Microwave Effects of Ethanol and Hexane-Mixed Solution as Revealed by In Situ Microwave Irradiation Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation.

Microwave heating is widely used to accelerate the organic synthesis reaction. However, the role of the nonthermal microwave effect in the chemical reaction has not yet been well characterized. The microwave heating processes of an ethanol-hexane mixed solution were investigated using in situ microwave irradiation nuclear magnetic resonance spectroscopy and molecular dynamics (MD) simulation. The temperature of the solution under microwave irradiation was estimated from the temperature dependence of the 1H chemical shifts (chemical shift calibrated (CSC)-temperature). The CSC-temperature increased to 58 °C for CH2 and CH3 protons, while it increased to 42 °C for OH protons during microwave irradiation. The CSC-temperature of CH2 and CH3 protons reflects the bulk temperature of solution by the thermal microwave effect. The lower CSC-temperature of the OH proton can be attributed to a nonthermal microwave effect. MD simulation revealed that electron dipole moments of OH groups ordered along the oscillated electric field decreased the entropy by absorbing microwave energy and simultaneously increased the entropy by dissipating energy to the solution as the thermal and nonthermal microwave effect. Ordered polar molecules interact to increase hydrogen bonds between OH groups as the nonthermal microwave effect, which explains the lower CSC-temperature of the OH protons. The nonthermal microwave effects contribute to the intrinsic acceleration of the organic reaction.

Aggregability of ß(1→4)-linked glucosaminoglucan originating from a sulfur-oxidizing bacterium Thiothrix nivea.

ß-1,4-glucosaminoglucan (GG) was prepared from the sheath of a sulfur-oxidizing bacterium Thiothrix nivea. Recently, GG was found to be adsorbed by cellulose (paper) and is therefore potentially applicable as an aminating agent for cellulose. We attempted to increase the yield of GG using a fed-batch cultivation method. Furthermore, the behavior of GG molecules in water was theoretically and experimentally investigated. NMR analysis in combination with molecular dynamics calculation suggested that GG molecules tend to form soluble aggregates in water. It was experimentally revealed that the self-aggregation is enhanced by the addition of NaCl and reduced temperature. Adsorption of GG onto cellulose via hydrogen bonding was confirmed by molecular dynamics simulation. Adsorption was also promoted in the presence of NaCl but was inhibited by a reduction in temperature. Only 11% of the amino groups in the GG-treated paper was reactive, suggesting that GG molecules adsorbed by the paper were forming aggregates.

Separation of D-amino acid-containing peptide phenylseptin using 3,3'-phenyl-1,1'-binaphthyl-18-crown-6-ether columns.

Several D-amino acid-containing peptides (DAACPs) with antimicrobial, cardio-excitatory, or neuronal activities have been found in several species. Here, we demonstrated the chiral separation of the antimicrobial peptide diastereomers, D-phenylseptin and L-phenylseptin using (S) and (R) 3,3'-phenyl-1,1'-binaphthyl-18-crown-6-ether columns (CR-I (+) and CR-I (-), respectively) and also investigated the underlying mechanism. First, using D-amino acid-containing tripeptide Phe-Phe-Phe-OH, we found that CR-I (+) could be used to recognize diastereomeric tripeptides containing an L-amino acid as the first residue. On the contrary, CR-I (-) enabled separation of a series of diastereomers with D-amino acid as the first residue. Therefore, we achieved separation of the stereoisomers using the chiral columns depending on the position of the D- amino acid in the peptide and demonstrated the orthogonality of separations of the chiral columns. Then, using CR-I (+), we separated amphibian antimicrobial peptide diastereomers, L- and D-phenylseptin, which have the sequences, L-Phe-L-Phe-L-Phe and L-Phe-D-Phe-L-Phe at their N-termini, respectively. In order to understand the host-guest interactions, we performed molecular dynamics simulations for L-Phe-L-Phe-L-Phe tripeptide-CR-I molecule complex systems. Three hydrogen bonds between the N-terminal amine group -NH3+ and the crown ether oxygens were the dominant interactions. The hydrophobic interactions between phenyl-rings in the chiral selector unit of CR-I (+) and the side chains of 2nd and 3rd residues of the peptide also contributed to the affinity. Our results show that the CR-I (+)-column can be applied for the separation of endogenous DAACPs generated by the post-translational modification.

DFT study of the influence of acetyl groups of cellulose acetate on its intrinsic birefringence and wavelength dependence.

The effect of the acetyl groups of cellulose acetate (CA) on its intrinsic birefringence and its wavelength dependence was investigated using density functional theory (DFT). Seven types of CA repeating-unit models that differ in their degree of substitution (DS) and substitution sites were used in the calculations. The results suggested that the intrinsic birefringence (Δn°) and its wavelength dependence significantly depended on the conformations of the acetyl group at C6. Additionally, the intrinsic birefringence of CA films was estimated as the ensemble average of the calculated Δn° values of the conformers. The increase in the DS of CA led to a more negative intrinsic birefringence and a larger wavelength dependence. The computational results were in good qualitative agreement with the experimental results and suggested that conformational variety and/or its control would be important factors for the design of optical films containing CA.

Retinal Configuration of ppR Intermediates Revealed by Photoirradiation Solid-State NMR and DFT.

Pharanois phoborhodopsin (ppR) from Natronomonas pharaonis is a transmembrane photoreceptor protein involved in negative phototaxis. Structural changes in ppR triggered by photoisomerization of the retinal chromophore are transmitted to its cognate transducer protein (pHtrII) through a cyclic photoreaction pathway involving several photointermediates. This pathway is called the photocycle. It is important to understand the detailed configurational changes of retinal during the photocycle. We previously observed one of the photointermediates (M-intermediates) by in situ photoirradiation solid-state NMR experiments. In this study, we further observed the 13C cross-polarization magic-angle-spinning NMR signals of late photointermediates such as O- and N'-intermediates by illumination with green light (520 nm). Under blue-light (365 nm) irradiation of the M-intermediates, 13C cross-polarization magic-angle-spinning NMR signals of 14- and 20-13C-labeled retinal in the O-intermediate appeared at 115.4 and 16.4 ppm and were assigned to the 13-trans, 15-syn configuration. The signals caused by the N'-intermediate appeared at 115.4 and 23.9 ppm and were assigned to the 13-cis configuration, and they were in an equilibrium state with the O-intermediate during thermal decay of the M-intermediates at -60°C. Thus, photoirradiation NMR studies revealed the photoreaction pathways from the M- to O-intermediates and the equilibrium state between the N'- and O-intermediate. Further, we evaluated the detailed retinal configurations in the O- and N'-intermediates by performing a density functional theory chemical shift calculation. The results showed that the N'-intermediate has a 63° twisted retinal state due to the 13-cis configuration. The retinal configurations of the O- and N'-intermediates were determined to be 13-trans, 15-syn, and 13-cis, respectively, based on the chemical shift values of [20-13C] and [14-13C] retinal obtained by photoirradiation solid-state NMR and density functional theory calculation.

Acetylacetonato-based pincer-type nickel(ii) complexes: synthesis and catalysis in cross-couplings of aryl chlorides with aryl Grignard reagents.

In this work, three different types of acetylacetonato-based pincer-type nickel(ii) complexes (2) were prepared. Complex 2a possessed the tridentate ONN ligand, which was constructed by the condensation reaction of acetylacetone with N,N-diethylethylenediamine. Complex 2b contained the PPh2 donor group in contrast to the NEt2 group in 2a, i.e., an ONP ligand framework. Complex 2c was composed of the NNN ligand, which was prepared by the reaction of 4-((2,4,6-trimethylphenyl)amino)pent-3-en-2-one with N,N-diethylethylenediamine. In addition to X-ray diffraction analysis, these complexes were characterized spectroscopically. Their catalytic activity for a cross-coupling reaction of aryl halides with aryl Grignard reagents was also evaluated. Among these complexes, 2b acted as an effective catalyst for the cross-coupling reaction using aryl chlorides as electrophiles. The electronic properties of these Ni(ii) complexes were investigated by cyclic voltammetry and density functional theory calculations.

The role of d-allo-isoleucine in the deposition of the anti-Leishmania peptide bombinin H4 as revealed by 31P solid-state NMR, VCD spectroscopy, and MD simulation.

Bombinin H4 is an antimicrobial peptide that was isolated from the toad Bombina variegata. Bombinin H family peptides are active against gram-positive, gram-negative bacteria, and fungi as well as the parasite Leishmania. Among them, bombinin H4 (H4), which contains d-allo-isoleucine (d-allo-Ile) as the second residue in its sequence, is the most active, and its l-isomer is bombinin H2 (H2). H4 has a significantly lower LC50 than H2 against Leishmania. However, the atomic-level mechanism of the membrane interaction and higher activity of H4 has not been clarified. In this work, we investigated the behavior of the conformations and interactions of H2 and H4 with the Leishmania membrane using 31P solid-state nuclear magnetic resonance (NMR), vibrational circular dichroism (VCD) spectroscopy, and molecular dynamics (MD) simulations. The generation of isotropic 31P NMR signals depending on the peptide concentration indicated the abilities of H2 and H4 to exert antimicrobial activity via membrane disruption. The VCD experiment and density functional theory calculation confirmed the different stability and conformations of the N-termini of H2 and H4. MD simulations revealed that the N-terminus of H4 is more stable than that of H2 in the membrane, in line with the VCD experiment data. VCD and MD analyses demonstrated that the first l-Ile and second d-allo-Ile of H4 tend to take a cis conformation. These residues function as an anchor and facilitate the easy winding of the helical conformation of H4 in the membrane. It may assist to quickly reach to the threshold concentration of H4 on the Leishmania membrane. This article is part of a Special Issue entitled: d-Amino acids: biology in the mirror, edited by Dr. Loredano Pollegioni, Dr. Jean-Pierre Mothet and Dr. Molla Gianluca.

Evaluation of hydrogen bond networks in cellulose Iß and II crystals using density functional theory and Car-Parrinello molecular dynamics.

Crystal models of cellulose Iß and II, which contain various hydrogen bonding (HB) networks, were analyzed using density functional theory and Car-Parrinello molecular dynamics (CPMD) simulations. From the CPMD trajectories, the power spectra of the velocity correlation functions of hydroxyl groups involved in hydrogen bonds were calculated. For the Iß allomorph, HB network A, which is dominant according to the neutron diffraction data, was stable, and the power spectrum represented the essential features of the experimental IR spectra. In contrast, network B, which is a minor structure, was unstable because its hydroxymethyl groups reoriented during the CPMD simulation, yielding a different crystal structure to that determined by experiments. For the II allomorph, a HB network A is proposed based on diffraction data, whereas molecular modeling identifies an alternative network B. Our simulations showed that the interaction energies of the cellulose II (B) model are slightly more favorable than model II(A). However, the evaluation of the free energy should be waited for the accurate determination from the energy point of view. For the IR calculation, cellulose II (B) model reproduces the spectra better than model II (A).

The role of CH/π interactions in the high affinity binding of streptavidin and biotin.

The streptavidin-biotin complex has an extraordinarily high affinity (Ka: 1015mol-1) and contains one of the strongest non-covalent interactions known. This strong interaction is widely used in biological tools, including for affinity tags, detection, and immobilization of proteins. Although hydrogen bond networks and hydrophobic interactions have been proposed to explain this high affinity, the reasons for it remain poorly understood. Inspired by the deceased affinity of biotin observed for point mutations of streptavidin at tryptophan residues, we hypothesized that a CH/π interaction may also contribute to the strong interaction between streptavidin and biotin. CH/π interactions were explored and analyzed at the biotin-binding site and at the interface of the subunits by the fragment molecular orbital method (FMO) and extended applications: PIEDA and FMO4. The results show that CH/π interactions are involved in the high affinity for biotin at the binding site of streptavidin. We further suggest that the involvement of CH/π interactions at the subunit interfaces and an extended CH/π network play more critical roles in determining the high affinity, rather than involvement at the binding site.

Application of the fragment molecular orbital method analysis to fragment-based drug discovery of BET (bromodomain and extra-terminal proteins) inhibitors.

The molecular interactions of inhibitors of bromodomains (BRDs) were investigated. BRDs are protein interaction modules that recognizing ε-N-acetyl-lysine (εAc-Lys) motifs found in histone tails and are promising protein-protein interaction (PPI) targets. First, we analyzed a peptide ligand containing εAc-Lys to evaluate native PPIs. We then analyzed tetrahydroquinazoline-6-yl-benzensulfonamide derivatives found by fragment-based drug design (FBDD) and examined their interactions with the protein compared with the peptide ligand in terms of the inter-fragment interaction energy. In addition, we analyzed benzodiazepine derivatives that are high-affinity ligands for BRDs and examined differences in the CH/π interactions of the amino acid residues. We further surveyed changes in the charges of the amino acid residues among individual ligands, performed pair interaction energy decomposition analysis and estimated the water profile within the ligand binding site. Thus, useful insights for drug design were provided. Through these analyses and considerations, we show that the FMO method is a useful drug design tool to evaluate the process of FBDD and to explore PPI inhibitors.

Dynamic Structure and Orientation of Melittin Bound to Acidic Lipid Bilayers, As Revealed by Solid-State NMR and Molecular Dynamics Simulation.

Melittin is a venom peptide that disrupts lipid bilayers at temperatures below the liquid-crystalline to gel phase transition temperature (Tc). Notably, the ability of melittin to disrupt acidic dimyristoylphosphatidylglycerol (DMPG) bilayers was weaker than its ability to disrupt neutral dimyristoylphosphatidylcholine bilayers. The structure and orientation of melittin bound to DMPG bilayers were revealed by analyzing the 13C chemical shift anisotropy of [1-13C]-labeled melittin obtained from solid-state 13C NMR spectra. 13C chemical shift anisotropy showed oscillatory shifts with the index number of residues. Analysis of the chemical shift oscillation properties indicated that melittin bound to a DMPG membrane adopts a bent α-helical structure with tilt angles for the N- and C-terminal helices of -32 and +30°, respectively. The transmembrane melittin in DMPG bilayers indicates that the peptide protrudes toward the C-terminal direction from the core region of the lipid bilayer to show a pseudotransmembrane bent α-helix. Molecular dynamics simulation was performed to characterize the structure and interaction of melittin with lipid molecules in DMPG bilayers. The simulation results indicate that basic amino acid residues in melittin interact strongly with lipid head groups to generate a pseudo-transmembrane alignment. The N-terminus is located within the lipid core region and disturbs the lower surface of the lipid bilayer.

Discrimination between naphthacene and triphenylene using cellulose tris(4-methylbenzoate) and cellulose tribenzoate: A computational study.

The mechanisms of naphthacene and triphenylene discrimination using commercially available cellulose tris(4-methylbenzoate) (CMB) and cellulose tribenzoate (CB) chiral stationary phases were investigated using molecular mechanics calculations. Naphthacene and triphenylene could be separated by liquid chromatography on CMB and CB, with triphenylene being eluted earlier than naphthacene on both phases. However, the corresponding separation factor is much larger for CMB than for CB. The docking of these polycyclic aromatic hydrocarbons to the above polymers suggested that the most important sites of CMB and CB for interacting with these hydrocarbons are located at equivalent positions, featuring a space surrounded by main chain glucose units and benzoyl side chains. The difference of hydrocarbon stabilization energies with CMB and CB agreed well with the observed chromatographic separation factors.

Achiral Molecular Recognition of Aromatic Position Isomers by Polysaccharide-Based CSPs in Relation to Chiral Recognition.

Chromatographic separation of several sets of aromatic position isomers on three cellulose- and one amylose-based chiral stationary phases was performed to evaluate the potential of a polysaccharide-based chiral stationary phase (CSP) in the separation of isomeric or closely similar molecules, and to understand the interaction mechanism of this type of CSP with analytes. Their ability of molecular recognition was quite outstanding, but the selection rule was particular to each polysaccharide derivative. In the series of analytes, cellulose tris(4-methylbenzoate) and tris(3,5-dimethylphenylcarbamate) exhibited a contrasting selection rule, and the recognition mechanism was considered based on the computer-simulation of the former polymer.

Ab initio studies on the structure of and atomic interactions in cellulose III(I) crystals.

The crystal structure of cellulose III(I)was analyzed using first-principles density functional theory (DFT). The geometry was optimized using variable-cell relaxation, as implemented in Quantum ESPRESSO. The Perdew-Burke-Ernzerhof (PBE) functional with a correction term for long-range van der Waals interactions (PBE-D) reproduced the experimental structure well. By using the optimized crystal structure, the interactions existed among the cellulose chains in the crystal were precisely investigated using the NBO analysis. The results showed that the weak bonding nature of CH/O and the hydrogen bonding occur among glucose molecules in the optimized crystal structure. To investigate the strength of interaction, dimeric and trimeric glucose units were extracted from the crystal, and analyzed using MP2 ab initio counterpoise methods with BSSE correction. The results estimated the strength of the interactions. That is, the packed chains along with a-axis interacts with weak bonding nature of CH/O and dispersion interactions by -7.50 kcal/mol, and two hydrogen bonds of O2HO2…O6 and O6HO6…O2 connect the neighboring packed chains with -11.9 kcal/mol. Moreover, FMO4 calculation was also applied to the optimized crystal structure to estimate the strength of the interactions. These methods can well estimate the interactions existed in the crystal structure of cellulose III(I).

Structure and orientation of antibiotic peptide alamethicin in phospholipid bilayers as revealed by chemical shift oscillation analysis of solid state nuclear magnetic resonance and molecular dynamics simulation.

The structure, topology and orientation of membrane-bound antibiotic alamethicin were studied using solid state nuclear magnetic resonance (NMR) spectroscopy. (13)C chemical shift interaction was observed in [1-(13)C]-labeled alamethicin. The isotropic chemical shift values indicated that alamethicin forms a helical structure in the entire region. The chemical shift anisotropy of the carbonyl carbon of isotopically labeled alamethicin was also analyzed with the assumption that alamethicin molecules rotate rapidly about the bilayer normal of the phospholipid bilayers. It is considered that the adjacent peptide planes form an angle of 100° or 120° when it forms α-helix or 310-helix, respectively. These properties lead to an oscillation of the chemical shift anisotropy with respect to the phase angle of the peptide plane. Anisotropic data were acquired for the 4 and 7 sites of the N- and C-termini, respectively. The results indicated that the helical axes for the N- and C-termini were tilted 17° and 32° to the bilayer normal, respectively. The chemical shift oscillation curves indicate that the N- and C-termini form the α-helix and 310-helix, respectively. The C-terminal 310-helix of alamethicin in the bilayer was experimentally observed and the unique bending structure of alamethicin was further confirmed by measuring the internuclear distances of [1-(13)C] and [(15)N] doubly-labeled alamethicin. Molecular dynamics simulation of alamethicin embedded into dimyristoyl phophatidylcholine (DMPC) bilayers indicates that the helical axes for α-helical N- and 310-helical C-termini are tilted 12° and 32° to the bilayer normal, respectively, which is in good agreement with the solid state NMR results.

Computational study to evaluate the birefringence of uniaxially oriented film of cellulose triacetate.

The intrinsic birefringence of a cellulose triacetate (CTA) film is evaluated using the polarizability of the monomer model of the CTA repeating unit, which is calculated using the density functional theory (DFT). Since the CTA monomer is known to have three rotational isomers, referred to as gg, gt, and tg, the intrinsic birefringence of these isomers is evaluated separately. The calculation indicates that the monomer CTA with gg and gt structures shows a negative intrinsic birefringence, whereas the monomer unit with a tg structure shows a positive intrinsic birefringence. By using these values, a model of the uniaxially elongated CTA film is constructed with a molecular dynamics simulation, and the orientation birefringence of the film model was evaluated. The result indicates that the film has negative orientation birefringence and that its value is in good agreement with experimental results.

Interaction of extracellular loop II of κ-opioid receptor (196-228) with opioid peptide dynorphin in membrane environments as revealed by solid state nuclear magnetic resonance, quartz crystal microbalance and molecular dynamics simulation.

κ-Opioid receptor is a member of the opioid receptor family and selectively interacts with the opioid peptide dynorphin. Extracellular loop II (ECL-II) of the κ-opioid receptor displays an amphiphilic helix in membrane environments and the N-terminal α-helix of dynorphin A(1-17) (hereafter DynA17) is inserted into the membrane with the tilt angle of 21° to the bilayer normal. ECL-II peptides (1-33), corresponding to 196-228 of κ-opioid receptor with [1-(13)C]- or [3-(13)C]-labeled amino acids were incorporated into large [dimyristoylphosphatidyl choline (DMPC)/ dihexanoylphosphatidyl choline (DHPC) = 3, q = 3] and small bicelle (q = 1) systems. (13)C direct detection with dipolar decoupling and magic angle spinning (DD-MAS) nuclear magnetic resonance (NMR) spectra were recorded, and the (13)C chemical shift perturbation clearly indicated that DynA17 interacts with ECL-II at the location of Val10-Ala15. Quartz crystal microbalance measurements were performed to determine the binding constant of ECL-II with DynA17 and indicated that the binding constant between DynA17 and ECL-II embedded in the lipid layer was 72 times larger than that between DynA17 and the lipid. The result of the molecular dynamics simulation clearly indicates that the C-terminus of DynA17 interact with the amino acid residues of the region between Val10-Gln14 of ECL-II. These results suggest that DynA17 interacts with the ECL-II of the κ-opioid receptor through a hydrophobic and short-lived electrostatic interaction with high affinity in the outer surface of the membrane.

Investigation of the structure and interaction of cellulose triacetate I crystal using ab initio calculations.

The crystal structure of cellulose triacetate I (CTA I) was investigated using first-principles density functional theory (DFT) calculations. The results are in good agreement with the experimental structure obtained by Sikorski et al. when performing the calculation with inclusion of the dispersion correction. However, the cell parameters calculated with inclusion of the dispersion correction are slightly smaller than those experimentally obtained, especially along the a-axis. This smaller cell parameter could be reasonably explained by considering thermal expansion effects, since optimization with the density functional calculation gives the structure without inclusion of thermal effects. The atoms-in-molecules (AIM) theory is also employed to identify and characterize interatomic interactions in the CTA I crystal. CH/O interactions sites are shown to exist in the crystal structure of CTA I. Moreover, CH/O interactions are considered the main interactions in operation to maintain the crystal structure of CTA I.

Computational study on the interactions and orientation of monoclonal human immunoglobulin G on a polystyrene surface.

Having a theoretical understanding of the orientation of immunoglobulin on an immobilized solid surface is important in biomedical pathogen-detecting systems and cellular analysis. Despite the stable adsorption of immunoglobulin on a polystyrene (PS) surface that has been applied in many kinds of immunoassays, there are many uncertainties in antibody-based clinical and biological experimental methods. To understand the binding mechanism and physicochemical interactions between immunoglobulin and the PS surface at the atomic level, we investigated the binding behavior and interactions of the monoclonal immunoglobulin G (IgG) on the PS surface using the computational method. In our docking simulation with the different arrangement of translational and rotational orientation of IgG onto the PS surface, three typical orientation patterns of the immunoglobulin G on the PS surface were found. We precisely analyzed these orientation patterns and clarified how the immunoglobulin G interacts with the PS surface at atomic scale in the beginning of the adsorption process. Major driving forces for the adsorption of IgG onto the PS surface come from serine (Ser), aspartic acid (Asp), and glutamic acid (Glu) residues.